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SYNOPSIS 

This thesis deals with dynamic stability of 
multi-span pipes conveying fluid where the flow velocity 
is either constant or a small harmonic component is super- 
imposed on it. The very general equation of motion is used 
to study the stability of the pipes. The element matrices 
are obtained by the finite element (Galerlcin) method. For 
steady flow stiffness, damping and mass matrices are 
obtained by finite element methods and a standard dynamical 
matrix is obtained. The critical velocities are obtained 
by solving eigen value problem. For harmonically perturbed 
flow the bounds for the principal primary regions of ins- 
tability are determined by Bolotin's method. It is shown 
that critical velocities and instability regions can be 
controlled by changing the position of the supports. 

Two general computer programs have been written 
to solve these problems. These programs are very flexible 
in nature and solution for various types of multi-span 
pipes can very easily be obtained by changing the input 
data . 



CHAPTER - X 


INTRODUCTION 


1.1 INTRODUCTION 

Pipes conveying fluids are encountered in various 
fields of engineering. To name some of them * oil pipe 
lines, propellent lines, pump discharge lines, heat 
exchanger tubes and coolant channels of nuclear reactors. 
An imperfect design of pipe lines may result in a leakage 
and unsatisfactory performance of the whole system. As 
the flow velocity is increased a certain velocity is 
reached when the pipe either buckles or flutters. With 
harmonic perturbations these systems are vulnerable to 
parametric instability also. 

In the recent years many pipe lines have been 
laid and this has accelerated the research work in this 
field. Lot of research has been done in the last thirty 
years. However most of literature is devoted to either 
single span pipe or periodic pipes. In actual practice 
however we come across nonperiodic multi -span pipes. 

The aim of the present work is to study the 
stability of multi'-span pipes conveying fluid at a 



constant velocity or when a small harmonic component is 
superimposed on it, by finite element method. All types 
of boundary conditions can easily be accommodated by this 
method. .The objective of present thesis is explained in 
detail in the last section of this chapter. 

2.2 REVIEW OF PREVIOUS WORK 

A detailed review of the literature with exten~ 
sive references is given in Paidoussis and Issid [2l] and 
Chen [6]. A brief review is given below. 

Interest in the subject of dynamics of pipes con- 
veying fluid was developed in 1950 in connection with study 
of vibration of Trans Arabian pipe by Ashley and Haviland 
[l]. They studied the problem of vibration of simply 
supported pipe. The same problem was studied by Housner 
[ 12 ] independently by usina a different approach. He 
found that at low fluid velocities the effect upon vibra- 
tions of the pipe was negligible but at certain high velo- 
cities pipe buckled like a column subjected to axial loading 

Long [l5] considered the case of cantilever pipe 
conveying fluid. He adopted an iterative procedure using 
a power series for mode shape, which was applicable to 
relatively small flow velocities and confirmed experi- 
mentally that forced motion of cantilever pipes was damped 
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by internal flow in the range of flow velocities considered. 

Benjamin [2] dealing with the general dynamical 
problem of cantilevered system of articulated pipes con- 
veying fluid was first to anticipate the phenomenon of 
instable oscillations (flutter) . He produced complete 
theory supported by experiments for articulated pipe sys- 
tems. He showed that when the system was vertical both 
oscillation and buckling instabilities are possible. 

However in general when motion is confined to horizontal 
plane, that is gravity is insignificant, buckling can not 
occur. Later Gregory and Paidoussis [lo] confirmed both 
theoretically and experimentally that at sufficiently high 
velocity cantilever pipes are subjected to oscillatory 
instabilities. Next Paidoussis C 1 9 ] included the effect 
of gravity and showed that in case of vertical cantilevers 
buckling instability is nob possible. This contradiction 
was cleared by Paidoussis and Deskins [£o].. 

Naguleswaran and William [17] studied 'both theo- 
retically and experimentally the effect of fluid pressure 
on the dynamics of the pipe. They showed that pipes with 
both ends supported may buckle at small flow velocities 
by the action of high internal pressure. 

Chen [5] studied the stability of pipes conveying 
fluid with upstream end fixed and downstream end supported 
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by displacement springs , so that boundary conditions are 
intermediate between clamped free and clamped~pinned . He 
showed that both buckling and oscillary instabilities are 
possible depending upon the spring constant. 

Paidoussis and Issid [ 2 l] derived the general 
equation of motion. They included the effect of axial 
movement of the pipe. They studied the dynamic stability 
of the pipe where the flow velocity is either constant 
or a small harmonic component is superimposed on it. They 
showed that with both ends fixed pipe is subjected not 
only to buckling instability but flutter instability is 
also possible at high velocities. 

Singh and Mallik [22] studied the wave propoga- 
tion and vibration response of a periodically supported 
pipe by propogation constant approach. A succeeding paper 
by same authors [23] was devoted to paiametric instabi- 
lities . 

Orris and Petyt [ 18 ] applied the finite element 
methods to propogation constant in periodic structures • 

Deb [7] applied finite element (Galerkin) method and 
obtained flutter and buckling instabilities of various 
types of single span pipes. 



5 


Kulkarni [ 13 ] studied instabilities of the periodic 
pipes due to parametric excitation considering the effect of 
axial movement of the pipe. He used the propogation cons- 
tant technique • 

1.3 OBJECTIVE AND SCOPE OP PRESENT WORK 

The aim of the present work is to develop finite 
element method for studying the dynamic stability of multi- 
span pipes conveying fluid. 

In second chapter element matrices are obtained 
by the finite element (Galerkin method ) for pipes convey- 
ing fluid when the flow velocity is either constant or a 
small harmonic component superimposed on it. For steady 
flow stiffness, damping and mass matrices are obtained by- 
the finite element method and a standard dynamical matrix 
obtained. Critical velocities are obtained by solving the 
eigen-value problem. For harmonically perturbed flow, 
bounds of principal primary instability regions are obtained 
by using Bolotin's method. 

In third chapter results for various types of 
pipe configurations are given and discussed. 

Conclusions are reported in the last chapter. 

The values or various matrices and listing of computer 
program are given in Appendixes . 



CHAPTER - II 


GENERAL FORMULATION 


This chapter deals with the finite element 
analysis of stability of pipes conveying fluid at cons 
tant velocity V and with a small harmonic component 
superimposed on it. 


2.1 EQUATION OF MOTION 


Consider a pipe conveying fluid at velocity V 
Fig. 1. Its equation of motion is, Paidoussis and 
Issid [2l] 

E*I — + El ------j + -j M^V 2 - T + PA (1 - 2vji) 

ax* t 9X 4 1 f 


where 

E* 

El 


[ (M. + m )g - M, S] (L - X) " 


3 W 

'f 1 “p /VJ “f 3t J ViJ " g 2 

+ 2 M f v me + <**, + V ^ ts + c n 


A 


+ (M^ + m ) 

£ P 3t 2 


= 0 


( 2 . 1 ) 


Coefficient of internal dissipation 
Flexural rigidity of pipe 
Mass of fluid per unit length 
Velocity of the fluid 


V 
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T Longitudinal tension 

P Fluid pressure above atmosphere 

v Poisson's ratio 

m Mass of pipe per unit length 

ir 

L Total length of the pipe 

c Coefficient of viscous damping 

ju,=0 For no axial constraint at the support 

ji.=l For axial constraint at the support 

• » • » (2.2) 


The equation (2.1) can be expressed in dimen- 
sionless form by defining the following quantities?. 


x = 


X 

L ' 


W = 


w 


, * = ( 


ei y 2 _i_ 

M f + m p l 2 


«.-C* 1 


M, 


1/2 jt M 1/2 & “f 

E(M f +m p )^ ^2 ' U = ( EI ) ' M f + m p 


y = 


"f * 1 ,.3, p „ Si! f „■ 

u g , l - E I ' r 


CL 


PAL 


EI 

2 


EI 


f>l(M f +m p )] 1/2 

(2.3) 


Substituting equation (2.3) into equation (2.1), the 
equation of motion becomes 


a 


,5 

<j w 

~— 

8x 5 T 




w w r z , _ 

+ ^ + J U - p + p (1 

8X 


2m) + (|3 1/2 & - 7)x 


2 

x(l- x ) j — ^ + 2(3 


, 1/2 


3 2 w 


3X 


u T*ri + y 


9w 

3X 


+ f 


T* 

. 9 2 w 

"H + TT2 


w 


o 


3T 
(2.4) 



8 


2 . 2 PIPES WITH STEADY FLOW 

For constant velocity, equation (2.4) becomes 


-P- + A + {u 2 -f* + p (1-2V«) - T 

3X 3p 3x * jpc 


+ 23 


1/2 


u 


3 2 v/ 
3x 3t 


+ 7 


' 5 V/ 
>x 


+ f 


3w 
3 T 


„ 2 
3 w 

2 

3T 


* 0 


(2.5 ) 


2.211 Finite Element Analysis: 

For finite element solution (Galerkin Method) 
of this one dimensional partial differential equation, 
the pipe is divided into number of finite elements. 
Figure 2, Typical finite element is shown in Figure 3. 
Its solution is assumed to be of the form of 
.(e) 


w 


( £ , T ) = L N a N, 


N 


J f w i 


<«2 


t 


or v/ e ^ ( £,t) = j^N^S) J £ w i ( tt ) J" 


(ne ) 


( 2 . 6 ) 


where 


Interpolating functions of the element 

w. Nodal parameters of the element 
3- / 

r Degrees of freedom of the element (to be 

"decided later) 

? Local co-ordinate of the element. 


Substituting equation (2.6) into equation (2.5), one 
gets the residue for the element, 
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( e ) 3 v/ 

R = a — ~ 

-r-4 


+ 7 (Xj + 5 ) 


-4 T (e ' 

ic 2 

+ 4— + 

f u 

3 5 4 

u 

3 2 w (e) „ 

3 v; ( 

95 2 

95 

3 V e) , 



• 8 2 w (e) 

P + p (l-2vn)-Y \ " •- -5 ■■ 

D €1^ 

+ 2P 1 / 2 u AilL 

3? a t 

(2,7) 


Minimizing this residue by the Galerkin method, 
rh R (e) d5 = 0 


or ^ Nh ( a 


9 5 w(e) 3 4 w (e) 

H 4 3T 


+ f u ~ P + p (l- 2 v/x ) 


, *2(e) n ( 

-? + 7 (x, + ;) i-V- + * — 

• r h 2 J H 2 35 


+ jfll/2 u iVfi 


■(e) 3 2 w ^ e ) 


) cIg, ' - o (2.8) 


959t DT 3t ^ 

In order to reduce the requirements of interpo- 
lating polynomials, the first four terms of equation 
(2.8) are integrated by parts, and one yets 


N . or. 


s 3 3 (e) h 3N. g 2 -(e) h h 3 2 N 

1 “i } - ~ a ^-j(ia— ) i + ; — ; 

3 E 3 A 35 sr 2 H n ar 


h h 31, -3 (e) 

1 1 r 1 „ a w j: 

{ 7‘ J 1 " ' Vjj -rr.H-il-L— H4,T-if _4_- 

O o 35 3«-,3 t 


+ N. 


-3 (e) h 3N, -2 (e) h , 3 N. -2 / i 

3 e 1 _ 1 3 w N I + 1 a w(e ) 


35 3 ^ 95 Q o 9^2 


a?- 1 95 


+ Kb .[ u - p + p (1-2V/J, )-7 


,(e ) h 


o y. 3 N . . ( e ) 

{U~ ~i> + p (1-2VJJ) - r>^ 2iL_ ae 


3w ' 

+ y x. N. -~- 

J 1 35 


(e ) h 


i- h 3N, -(e) 

I -TX, d < 


,(e) h 


J" — ^ d £ + rp N.- 

J O 95 ’ 1 


Contd 
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- ' i 


h 3N i 


a w 


(e) 




1/5 ft 2 ( e ) 

l 1 ' 2 u — — d£ 

1 dt D r 


N. 2/3- 


+ £ Nj 3W 


(e ) 


,2 (e) 
d w 


as + r* n. — as = o 

3t o 




(2.9) 


,(e ) 


Highest derivative of w with respect to < { in any 
integral in equation (2.9) is second. Thus interpolat- 

( g \ 

ing function N. should have compatibility of w and 
Highest derivative in equation (2.9) is third. 


9w 


(e ) 


3? 


Thus interpolating function N. should have completeness 

j3 w ( e ) 

and s; — . The function satisfy- 


of w 


(e) 9w (e) 3 2 w 8-V 


9? ' D? 2 


3 %r 


ing these requirements of compatibility and completeness 
is 

w (e) = a + bx + cx 2 4- dx 3 


or w 


( 0 ) 


= j (3-2^)^^^ , Z 2 2 (3-2 j''Wjj 


( 2 . 10 ) 


or w (e) * Ln (?) J {w ( t )} ne 
where = 1 - ?/h, $ 2 = C/h 

that is the typical element has four degrees of freedom. 

Using equation (2.10), the equation (2.9) can be 
written in the following matrix form 
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[m] (e) (w | (ne) + [c] (e) fw}- (ne) + [k] (e) fw} (ne) 


i(ne ) 


<! n a 


- 3 w 


- < N 


3 3 w (e) 

9£ 3 


(e ) h 
•) i 


9T 




3N a £ 


9w 


(e) 


*5 




9T7 


I 

o 




DM » 2 w (e) * 


3? 


35' 


2 C G ) ll 

~*jN(u -r+ P (1-2\*1) - r) J 

° J 


1 


r (e ) ‘ 

r N (x . + F , ) 3w 


J 


35 


I h 

r * 

J ° 


where 


h 

I 

o 


( 2 . 11 ) 


[m] (e) = / h { N } (_ N j dC (2.12) 

[c]^ = a / h {n"}{ N"Jd£ + 2I3 1 / 2 u 

+ f t { N / t N j (2.13) 

[k] (e) - 4 h fN"}.[ N"Jd5 

- (u 2 -p + p( 1-2'»M ) - r) <«$ h { N' J-[n' | CU 

“ rx j ^ i N '} |_ N ' J dK - r < h { N 1t N 'J d? (2.14) 

Values of matrices [m] (e) , [c] (e) and [k] -are 
given in Appendix I. Using the values of Ik from equa- 
tion (2.10) in the right side of equation (2.11), 


gets 


one 
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[m] (e) {w} <ne) + [ c ] (e> {w} (ne) + [k] (e) {w> (ne) 
‘ . 3 4 ,„< e > 


f 

r * 


~ a 


+4 


3? DT ■* j 
0 

^4 ( e ) 

3 w I 

9 5 3 3 X ' k 
0 


- a 


0 


.3 (e) 

3 V7 


»S 2 ar!j J 

1+ : 


o 


a. 


a 3 w (e) j 

35 2 9T k 


9 3 w 1 
H 3 I J 


o 


aV e) 


ar 


o 


o 

9 2 w (e) 


a? J 

o 



(u 2 - r +p(1-2vr)-7) ~ 

(e ) , ‘ 




95 Ij 


35 

j 


0 



0 


*+ * 

-~(u 2 -r +p(l-2vju)-T ) ~~ 


U-- 

Sl .(e) 

~rx— 

*7 

# 


a 5 

k 


H 

k . 


0 



0 


-i 


.w 





*2 < e > 
q W 

95 2 

J u ' ' j L 

(2.15) 

The equation (2.15J is for one finite element. These are 
assembled for the whole domain, 

[M] {w} (n) + [C] {w} tn > + [k] (wt- (n > = 
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After applying boundary conditions, the equation (2.16) 
reduces to 
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V 


Cm] 


BC 


•T w j 
I. J 


l (n) 


[c] 


BC 




1 (n) 


[K] 


BC 


/A (n) 

t w r 


For boundary conditions, see Appendix 2. 
2.2.2 Method of Solution: 


0 

(2.17 ) 


Equations (2.17) are set of homogenous differen- 
tial equations. The solution of this homogenous set of 
differential equations is obtained as detailed in 
Meirovitch ( 1 4 ] - But it may be noted that the matrices 
[m] 0 _, [c]__‘ and [k] need not be symmetric as stated 

in Meirovitch [l4], see Frazer, Duncan and Collar [ 9 ]. 
Using generalised velocity w as auxiliary variables, 
n second order ordinary differential equations (2.17) 
are converted to a set of 2n first order ordinary diffe- 
rential equations. 


[m] { )} + [K]|n(t)^ = J o] 


where 


< 1) (T.)>- 


and 


Cm] 


j -j W ^ 

j~ w (.t 


[ 0 ] [m] 


BC 


^bc^-^bc 


and [l<] 


-Cm] 

Co] 


BC 


(2.18) 


(2.19) 


Co] 

^ K ^BC 


( 2 . 20 ) 
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Now we have the differential equations which in 
the matrix form are 

[m] + [k] p = { oj (2.21) 

set •» T) (t; ) : - = e^ i 0 l (2.22) 

*’ J L i 

where 0 consists of a vector consisting of 2n constant 
terms. Substitution of equation (2.22) in equation (2.21) 
leads to the eigen-value problem 

¥ [m] •{ 0 > + [k] 0 } * 0 (2.23) 

which can be written as 

[d] { 0 j- = | { 0 .} (2.24) 

where ** 

r to] [i] 

(d] = - [k]- 1 [m] = _[k]- 1 Cc] <2 - 25) 

The eigen value of equation (2.23) will in general 
be complex. The real part of eigen value determines the 
stability of pipe and imaginary part gives the frequency 
of pipe. If the real part is positive and imaginary part 
is nonzero then there will be flutter. If real part is 
positive and imaginary part is zero the pipe will buckle. 
The velocity u at which pipe fails either due to buckl- 
ing or flutter is known as the critical velocity, u c> 

A library subroutine EIGRF was used to calculate 
the complex eigen values. 
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2.3 PIPES WITH HARMONICALLY PERTURBED FLOW 


Equation (2.4) is general governing differential 
equation of the pipe conveying fluid. 


Now a case is considered where the flow velocity 
is perturbed harmonically, that is 


u = u (1 + 6 cos Q x ) (2.26) 

where u q is mean flow velocity, 6 is excitation para- 
meter and Q is non-dimensional frequency given by — . 

y 

Substituting equation (2.26) in equation (2.4) 

one gets 


3 4- 

or. + [ u 2 (1+6 cos Qt -i ! + p (l-2vju) 

ax 9 t ax 


>r (i - x) - /3 1/2 u 


o 


.2 

6 Q (1 - x) sin R t] — 

9X 


+2P V 2uo ( i+« COS0T, ^C + ,g + f || 

,2 

+ = 0 (2,27) 

anr 


The coefficient of some derivatives are time 
varying, hence the pipe is parametrically excited. Because 
of parametric excitation, instability will occur over a 
range of frequencies. 


In order to obtain the primary instability 
regions,, we can write, 3olotin [ 4 ] 
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w ( x, T ) -- 2 (Y^ (x) sin (^ 2 ~) + 2^ (x) cos (^?p) 

k = 1,3,5 (2.28) 

The regions of primary instabilities can be 
obtained quite accurately by taking K = 1 approximation, 
Paidoussis and Issid [2l], Hence one can write 


w ( x , % ) = Y (x) sin + Z (x) cos (~ 4=p- ) (2.29) 

Substituting equation (2.29) in equation (2.27) 
and equating the coefficients of sin (^~) and cos 
separately to zero, one gets 

34 + ( u* u - 6+ -S) _p+ p (i 

dx 


•Q T 


\ c2 Y + * az 


» v jB 1 / 2 u <5 q 
3Y o [ (1 . x) 


a Q d Z f 0 

___ 

dx 


d z 


Z = 0 


dx 


4 


i 0 2 z + i 52 . pl/ S> 6ii 

? 0 Z + 1 35 2 


2VJ.I ) ~ 

r (l-x))A 

dx 

,2„ 
d Z 

dx 

(1 ~ 5‘ ) d£ J 


( 2 . 30a ) 

- 2vju) 

, 2_ 

~ r (1 - x))2-| 
dx^ 

? 

x d y 

9 

dx 

_ <i + h ftl] 

Kl h dx J 


+ ft ^ ft . Y , f S3 y — o 

dx 

Let us take 


( 2 . 30b ) 


~ : 1 n 2 
A 2 “ 4 0 


£ 

2 


A., = uf (1 — 6 + -%) - P + p (1 ■- 2v ) - y 


(2.31) 

(2.32) 
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A 3 


/3 1//2 u Q 


0 P 1//2 U &Q 

A -(1.2, 0__ 

A 4 ( 6 2 


(2.33) 

(2.34) 


A 5 = u 2 (1 + 6+ |~) - + p (1 - 2vju) - 7 


(1 + ~) 


„ P 1//2 u 6 Q 

^ \ o 


( 2 . 35 ) 


(2.36) 


c ' 2 

Using equations (2.31) to (2.36), equations (2. 


become 


+ A 1 ^ - A 2 y + 7 (X g) - A (1 - x) 
dx A 1 dx cix cix 3 dx ‘ 


+ A 


dZ a_fi d_Z £ Q 

4 3x " 2 . 4 ~ 2 

dx 


Z = 0 


(2.37a) 


+ A 5 ~ A 2 Z + T ix (x H ) - A 3 (1 


- x) 


l 8 * 


dx' 


+ A §1 + SLR ±1 + LJ1 Y - o 

+ A 6 dx + 2 . 4. 2 * U 

dx 


(2.37b) 


2.3.1 Finite Element Analysis s 

Equations (2.37) are coupled differential equa- 
tions. For finite element solution of these two coupled 
differential equations, the pipe is divided into number 
of finite elements, Figure 2. Solution over the finite 
element is assumed of the form of 

Y (e) U) = | N i (?) | |Y,J (ne) (2.38a) 


30) 
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and 

Z (e) ( 5 ) = j B. U) | j Z. L (ne) (2.38b) 

1 i 1 j 

where , 

Ni , Interpolating functions 

Nodal parameters 
S Local coordinate 


Substituting equation (2.38) in (2.37) one gets the 
residues and R^ e ^ for equation (2.37a) and (2.37b) 

respectively. 


R 


(e) d 4 Y (e) , . d 2 Y (e) 


d ? + Al d?" 


( e ) 

- a 0 Y v ; + r 


d_ 

d 0 


( (x . + 0) 


dY 


(e) 


d 0 


-A (1 - x » d 2 z (e> dZ (e > o_£ dfz f_H „(e) 

A 3 *j E) “TP" A 4 d 5 - 2 df 4 - 2 2 


• I « t 


<2. 39a) 


-.4(e) 4.2(e) , •. , ... „( e) 

R b ■ tnr- + A s hr- ~ A a 2 + T ((x j +E) ir 1 


dr 


dO 


a m v rl d 2 Y (e) . dY (e) a q d 4 Y (e) f Q, Y (e) . 

“ A -> (i - + A e di — - + ~r + ~T" Y 


do' 


(2.39b) 


.th 


where is the global coordinate of the j node 


Minimizing the residues by the Galerlcin method 
^ N i R a ^ = 0 (2.40a) 

i h B i R b^ S) d5 = 0 (2.40b) 
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Treatment for (2.40a) and (2.40to) is similar. 

So details are given only for (2.40a) and results of both 
are given in equation (2.46). 

( e ) 

Substituting the value of R from equation 

8. 

(2.39a) in the equation (2.40a), one gets 


i h + A ) 

d 5 


- A 3 (1 - x. - K) 


« o 


dV e) 

- A, Y (e) 


(U.+ 


d^ 2 

2 

d? 

J 

d? 

d 2 Z (e) , 

A < 3z(e ’ 

a Q d 4 

z (s) 

f £3 

d5 

A 4 3T 

dE 

4 

2 


* • 

• to 

(2 

.41) 


(e) 


Z ve ')d? 


In order to reduce the requirements of the enterpolating 
polynomials first, second, fourth, fifth and seventh terms 
are integrated by parts, and one gets 

.h ^i d 2 Y (e) . * ,h dN i dY (e) fh , 

l df ~2“ d5 “ A 1 & d? d? dC " A 2 i N i Y d5 


d?' 

,h dN i dY^ e 5 


h dN i dY (e ' 


•^jTd f3T d5 - r 4 E ac JT d? 


h dN i dz (e) 


+ A 3 (1 “ X j ) £ d ? 

, (e ) 


A 3 f K.. 


dZ 


- a 3 £ 

(e) 


. dw. dZ 1 '®^ 

h i 

d? d? 


d? 


, * r h dZ 

• air dr + A. / N. ■*■= — d? 
l d^ ^ 4 o id? * 


a 0 .h d N i d 2 Z (e) ., 

t~ 6 ~~~y dc 


dr 


d V 


N x Z (e) d? 


N . 


d 


3 Z (e) h dN i dV e) 


d r 


O 


de d?" 


+ A. N. 3y 
1 i d? 


dY^ e ^ h 
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- Xj - 5) q - 

z <e> * 

d * i 

- r h, 

i 

(x. + £ 

dY (e) 17 
’ d * ' o 

7 ( e ) h 

o!ll ! + 

d F, 6 

a Q df 'i 

~T~ dc 

d 2 z (e) 

d£ 2 

h 

' = o 

( 

O 

(2.42a) 


nr. Q 

-~5^ N. 

2 i 


Similar equation. (2.42b) is obtained from equa- 
tion (2.40b), Highest order of derivative in the inte- 

(e ) 

grals in equations (2.42) is two. Therefore Y and 

Z*' 0 ^ should have compatibility of (Y, and (Z, ^|) 

respectively. Highest order of derivatives in equations 

(2.42) is three. Therefore Y^ e ^ and Z ^ e ^ should have 

2 3 , 2d 

completeness of (Y, ™ £ , ^-—t?) and (Z, S , — r) 

dir as 3 as 2 as 3 

respectively. Thus 1st and can be chosen identically. 
Interpolating functions satisfying these requirements 
are. Able and Desai [s] 


N i U) j*l B i (C) J= \A (3 - *2 , (3 “ 2 h } r h hh 


vjhere 


E j = 1 - ^ and E 2 - jf 


(2.43) 


Using equation (2.43), the equation (2.42a) 


becomes 


<5 h {N" } |N" J d e { Y } (ne) - Aj /‘{N' } (_N'j dE (Y } (ne) 

- *2 ?{« }tM J« {Y} (ne> - .Xj J *{Y]f ns) 

- 7 A{ N ' } Id' J «{*]’ ne) + A 3 d - q> i h {N'}[N'Jd6[2j. (ne) 
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(2.44a) 
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Similarly equation (2,42b) can be written 


as 


f { N \}.L N "'J « { 2 } (ne) - A 5 g h {N*}iN'J <35 {Z} (ne) 

“ A 2 { N lt N > d S {2 ji ne} - r Xj ; h { N' }j_N' J d5 {z } (ne) 

^ S't N L N' J d 5 •• Z ne ^ + A,( 1 -x j ) ^ * N' ji- j N' fd£ *Y j ^ ne ^ 


- r 


- A 3 ,£ h 5{ N'} ^ H' l d{ JY J 


<ne) - a 3 g h (w >|n'J <ae *! y ne ! 


+ A 6 f5 h { S }L N ' j d « {Y} (ne) tijS^ {t»"}|N"JdU Y} <ne) 


r h 

f Q ,h 


+ & 11 i N L N J d? [y } ^ ne) 


f a 3 z (o) I 



d 5 3 

j 


0 


>+ < 

. d 3 z (e 



CO 

Its 

V 

! 

k 


0 



U'** 

L 


O 


d 2 Z (e) 


d£‘ 


O 


dV 6) i 

d£ f 


kJ 


:> + A s 5 


dZ 


(e ) 


d 5 


dZ (e > 

“ ~dT~l 


+ A 


- (1-x) 


(e ) 


d? 


■O 


- x 


d 5 


I 0 

L 

7 (e), 

I 


’j 


0 


3 5 


+ ( 1 ~x ) 


• (e) 


? + * 5 


d 5 


0 


+ x 


z 


(e) 


O 
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(2.44b) 
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Matrices [a], [b], [e] and [g] are given in 
Appendix 3 . 


The equations (2,45) are for one element. Now 
assembling it for whole domain and applying boundary 
conditions/ one gets 


or 



Various types of boundary conditions are listed 


in Appendix 4 
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2.3.2 Method of Solution: 

Bolotin has shown that instability regions can 
be obtained by equating determinant of matrix [p] to zero 
Hence frequencies 0 } , }. 1 , 2 are the frequencies 

where Det [f] = o and they give boundaries of instability 
regions . 

A Library subroutine PO 3AAF/NAG was used to cal- 
culate the determinent of matrix [f]. The instabilities 
regions are plotted in £3 - 6 


space » 



CHAPTER - III 


RESULTS AND DISCUSSION 

This chapter deals with numerical results obtained 
for pipe conveying fluid for various types of pipe confi- 
gurations. In Section 3.1 stabilities of pipes are dis- 
cussed when the fluid is flowing at constant velocity. 

In Section 3.2 instability regions for various types of 
pipe configurations are plotted when a small harmonic 
component is superimposed over the mean velocity. 

3.1 PIPES WITH STEADY FLOW 

In this section critical velocities of various 
types of pipe configurations have been found by the finite 
element methods. 

3.1.1 Number of Finite Elements 

To begin with one must decide the number of finite 
elements to be used in the analysis. For every type of 
boundary condition the frequencies were calculated for 
few velocities with increasing number of finite elements. 
We selected that number of elements beyond which the 
increase in number of finite element did not change the 
value of frequency significantly. 



29 


Table 1 shows -the frequencies of two equal span 
simply supported pipe with increasing number of finite 
elements for 0 = 0.5, u = 2 and a = = p = 7 = f =* 0. It 

is observed that the value does not change significantly 
after eight elements. So the number of finite elements 
taken for this case was eight. 

3.1.2 Results and Discussions 2 

Using the finite element equations developed in 
Chapter-II the critical velocities for various configura- 
tion of pipes are obtained. 

3. 1.2.1 One Span Pipes 

Complex frequencies for one span simply supported 
pipes were obtained by Paidoussis and Issid [2l] for 
/3 = 0.5, a=p~P = y ~ f = 0. The same problem was 
studied by Deb [ 7 ] by using finite element method. 
However Deb [ 7 ] did not study the effect of internal 
dissipation (a) and gravity ( 7 ) effects. 

Figure 4 shows effect of coefficient of internal 
dissipation on the same system. The value of a is 
taken equal to 0.005. It is found that critical 
velocity does not change but the symmetry about 
imaginary axis is lost. These results match with 
those obtained by Paidoussis and Issid 
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critical velocities of all these four cases are less 
than the critical velocity of three equal span pipe. 

3. 1.2. 3 Multi-span Pipes, One End Pi 2 d and Other 
Supports , Simple Supports 

Critical velocity for two equal span pipe, case (1) 
Table 3, for B = 0.5, a =f = p = y = f = 0, is 7.12. 
This critical velocity is more than the critical 
velocity of pipe when all the supports are simple 
supports as expected. Next this pipe is considered 
with two unequal span lengths, cases {2, 3), Table 3. 

It is found when the Intermediate support is shifted 
away from the fixed end critical velocity becomes 
7.42 which is more than that of the pipe considered 
in case (1). On the other hand if the intermediate 
support is shifted towards the fixed end the critical 
velocity becomes 5 .56 which is less than that of 
case (1 ) . 

Next three span pipes are considered, case (4, 

5, 6,7), Table 3. It is found for equal span pipe 
the critical velocity is 9,93. Critical velocities 
for unequal spans, cases (5, 6, 7) are 8.60, 8.10 
and 10.06 respectively. It may be noted that these 
critical velocities are more than those if all the 
supports are simple supports, see Table 2, 
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3. 1.2. 4 Multi-span Pipes with One End Fixed, Other 
End Free and Intermediate Supports Simple 
Supports 

Critical velocity for this two equal span pipe, 
case (1), Table 4 is 4.64 for P -0.5, cc =P = p = y 
= f = 0. It the position of intermediate support is 
shifted towards the free end the critical velocity 
becomes 5.65, case (2), Table On the other hand 
if intermediate support is shifted toward the fixed 
end, the critical velocity becomes 3.96 which is less 
than that of pipe considered in case (1). 

Next three • span pipes are considered, cases 
(4, 5, 6, 7) Table 4. The critical velocities for 
these cases are 6.48, 7.41, 5,79 and 5.24 respective ly . 
It may be noted that these pipes are similar to these 
considered in cases (4, 5, 6, 7) of Table 3, except 
that the right hand support has been removed and is 
free now. It is observed that if the right hand 
support is free, the critical velocities decrease. 

3. 1.2. 5 Multi-span Pipes with Both Ends, Fixed and 

Supports Simple Supports 

Critical velocity for two equal span pipe, case 
(1), Table 5, for 0 =0.5, a * [» a p = ‘Y = f = 0 , is 
9.0. For unequal span lengths case (2), Table 5 the 
critical velocity decreases to 8.65. 
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For three span pipes , cases (3, 4, 5, 6) critical 
velocities are 11.76, 9.90, 11,24 and 10.26 respec- 
tively. It should be noted that these pipes are 
similar to those considered in cases (4, 5, 6, 7) 

Table 3, except that the right hand support has been 
fixed. It is observed that in the present cases 
critical velocities increase. 

3. 1.2. 6 Multi-span pipes with One End Fixed 
and Other Supports Elastic Supports 

Table 6 gives the results for the pipes whose one 

end is fixed and all other supports are elastic 

supports. In cases (1, 2, 3, 4), Table 6 elastic 

supports are displacement springs. For 0 = 0.5, 

a=P=p = V = £ - 0 and spring constant a £ = 10.0 

the critical velocities are 8.5 2, 7.14, 7.08 and 

6.93 respectively. These velocities are less than 

that of those considered in cases (4, 5, 6, 7), 

Table 3. Where all the supports are simple supports 

as expected. 

Next pipes with torsional springs are considered, 
see cases (5, 6,7, 8), Table 6. The critical velo- 
cities are 13.41, 10.92, 10,71 and 10,38 respec- 
tively. These velocities are more than that of these 
considered in cases (4, 5, 6, 7), Table 3 as expected. 
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3. 1.2, 7 Multi-span Vertical Pipes 

Table 7 gives the critical velocities for vertical 
multi-span pipes. It is observed that critical velo- 
cities are more in case of downward flow than those 
for horizontal flow, and critical velocities are less 
in case of upward flow. It is observed that effect 
of gravity is more pronounced in multi-span pipes 
if one of the ends is free , 

3.2 PIPES WITH HARMONICALLY PERTURBED FLOW 

In this section we will discuss the stability of 
pipes when a small harmonic component is superimposed over 
the mean velocity. 

3.2.1 Results and Discussions s 

Using the finite element equations developed in 
Chapter-II the regions of principal primary instability 
are obtained for various types of pipe configurations. 

These regions are plotted in Q - 6 space. 

3. 2. 1.1 One Span Pipes 

In Figure 5 instability regions for one span 
simply supported pipe are plotted for u = 0.6JT, 
a =r* = p = y = f = 0 and = 0.4 and *=0.8 

These results match with those obtained in [21,13]. 
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1/2 

For 3 =0,8 results are also given in Table 8 and 

compared with Kulkarni [ 1 3 ] * Excellent matching is 
observed, it is observed if the value of 3 is increased 
instability region broadens. 

Results for - fixed-fixed pipe are also given in 

the Figure 5. Instability regions are plotted for 
1/2 

P = 0.4, a =p =p= y = f = 0, u = 0.6 rr and n. It 
is observed that with increase in velocity u the 
instability region shifts to lower value of Q. 

In Figure 6 instability regions are plotted for 
fixed simply supported pipe for 3^ // ^ = 0.4, a = 

= p ~ y - f = Q, u = 0.6 tt and it . Again it was 
observed that with increase in value of u instabi- 
lity region shifts to lower value of Q. In the same 
figure instability regions for cantilever pipe are 
also plotted for ~ 0.4, a = P =3= r = f = 0, 

u = tt and u = 1.2 tt. A striking feature of this 
instability region is that it starts after a finite 
value of 6 . Again instability region shifts itself 
to lower value of Q for higher velocity. 

3. 2. 1.2 Multi-span Simply Supported Pipes 

Instability regions associated with first two 
modes of two equal span simply supported pipe for 
u = 1.2 TT, 3 1//2 = 0 .8, a - p = r - f ~ 0 are 
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given in Figure 7 and Table 9. The results obtained 
by Kulkarni [ 1 3 ] are also given in the same • figure 
and table. By present analysis we obtained broader 
instability regions. The difference arises due to 
the change in the governing differential equation of 
motion, Kulkarni [13] has taken L in third term in 
equation (2.1) as length of the sy m, where actually 

L is the ovex'all length of the pipe. 

Figure 8 and Figure 9 give the instability regions 

associated with first two inodes -of two and three span 

1/2 

(equal and unequal) simply supported pipe for £ ' =0.4/ 
or. — r* = p = r = f = 0, respectively. For two span 
pipe velocity u is 0.6 v and for three span pipe it 
is 1,8 it. It is seen that for unequal span lengths 
the instability region associated with first mode 
shifts to lower value of Q in comparison to first 
mode of equal span pipes . On the other hand the 
instability region corresponding to second mode shifts 
to higher value of Q in comparison to instability 
region corresponding to second mode of equal span 
pipes. The instability region corresponding to first 
mode of one span pipe is also plotted in Figure 8. 

It can be seen that for one span pipe the instability 
region are at very low value of Q, 
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3.2.1 .3 Multi-span Pipes, one End Fixed and 
Other supports Simple Supports 

Figure 10 and Figure 11 show the instability 
regions associated with first two modes of two and 
three span (equal and unequal) pipe for = 0.4, 

a. =p = p s= r = f = o respectively. For two span 
pipe u is 1.2 it and for three span pipe it is 1.8 ft. 

It is observed that regions of instability associated 
with first mode of unequal spans are at lower value 
of Q in comparison to that of equal spans. On the 
other hand the instability regions associated with 
second mode are at higher value of 0 in case of un- 
equal spans. In Figure io instability regions corres- 
ponding to first mode of one span fixed - simply 
supported are also plotted. It is observed that for 
one span pipe instability region is at very low value 
of Q . 

3. 2. 1.4 Multi-span Pipes, Both Ends Fixed and 
Other Supports Simple Supports 

Figure 12 gives instability regions associated 
with first two modes of a two span (equal and unequal) 
pipe for u = 1.2 n, 8 1 / 2 = 0.4, a = p = p- 7' = f = 0. 
The region associated with first mode is at lower 
value of Q for unequal span length. But region of 
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instability associated with second mode is at higher 
value of Q for unequal spans. The results for one 
span are also plotted for comparison purposes . 

Figure 13 gives instability regions associated 
with first two modes of three span (equal and unequal) 
pipe for u » 1.8 h, 0.4/ ot =f*=p = y » f = 0. 

Here the instability regions associated with both the 
modes are at lower value of Q for unequal spans. 

3. 2. 1.5 Multi-span Pipes, One End Fixed, Other and 

Free and Intermediate Supports Simple Supports 

Figure 14 shows the instability regions for a 
two span (equal and unequal) pipe. The interesting 
feature is that the instability region starts from 
a finite value of 6 and both modes are at higher 
value of Q for unequal spans . 

Figure 15 gives instability regions for a three 
span pipe for u = 1.8 n, /3 1//2 = 0.4, a si" 1 =p=r=f=0. 

The .instability regions of both modes are at lower 
value of 0 for unequal spans. The instability region 
corresponding to the second mode brodens in case of 
unequal spans , 



CHAPTER IV 


CONCLUSIONS 

Based on the results obtained in Chapter III follow- 
ing conclusions are drawn. 

1. The results by finite element method for multi- 
equispan simply supported pipes with steady flow match with 
those obtained by earlier authors. This shows that finite 
element equations of motion are correct and can be applied 
for multi-span (equal and unequal) pipes with any type of 
supports . 

2. The results by the finite elements method for multi- 
equispan simply supported pipes with harmonically perturbed 
flow match with those obtained by earlier authors for ( 3 = 0 . 
For (3^0 results differ from earlier results because of 
correct equation of motion used here. Thus the finite ele- 
ment equations of motion developed here are correct aed can 
be applied for multi-span (equal and unequal) pipes with any 
type of supports . 

3. The method developed is very general and flexible. 
The critical velocities and instability regions for any type 
of pipe configurations can be obtained by changing the input 
data in the computer program. 
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4. The column matrices, N(u z - p +p(l-2vjx) ~ 7 ) ^ I i 1 

1 q 

and') v N(x. +H)^~ !■ of equation (2.11) should be trans- 

J ^ 6 •' 

f erred to square matrix of left hand side, similarly, third, 
fourth fifth column matrics on right hand side of equation 
(2.45) should be transferred to square matrices of left hand 
sides, as suggested by Deb [ 7 ]. These column matrices can 
not be neglected as in suggested by the work of Szabo and 
Lee [ 24 ]. 

5. Critical velocities of pipes can be controlled by 
changing the position of the supports. Critical velocities 
of vertical pipes decrease for upward flow and increase for 
downward flow in comparison to horizontal flow. This 
effect is more pronounced if one end is free. 

6. The principal primary instability regions broadens 
when the value of 8 is increased. The principal primary 
instability regions shifts to lower value of Q for increase 
In velocity u. 

7. If one end of the pipe is free the principal 
primary instability regions starts from a finite value of 

6 '*• 

8. Instability regions can be controlled by changing 
the position of the supports . 
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TABLE 2 
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TABLE 3 


Critical velocities of Multi-span pipes - o ne end f i xed 
and other supports si mple support 

3 = 0.5, a = p = p = r = f = o 


Type of pipe 


4- 1 


r 


~r~ 

2 * ■ ■" 
“I 


•L/2 




-L/2 


"h 


Xi/4 — 


-3L/4- 


~ZT 


3 


3 1 


— 3L/4- 




j*~L/ 4~*j 
^ 


;Xi*~“1j/ 3 — — -uMi- 1>/ 3— ™<im«*—L/ 3 • " t»j 

4 zr "S 



E 

u* 


~L/ 3 — —At/yii/ 6 J pf — L/ 2 — — 

zk zr~' 

~ j,/ 2 Ii/ 3 — — 

~ ~ 7T~ ~ZS 'A 


Number of ! 

Element i 
1 
I 

Critical 

Velocity 

8 ! 

J 

7.12 

8 

.5 *56 

8 

7.42 

> 

12 

9.93 

12 

9.60 

12 

8.10 

12 

10.06 

• M 
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f° r M ulti-span pip es with_both end fix ed 
and o t her support s sim ple su ppo rts ’ 

£ = 0.5, a = p *p=ysf = o 


Type of pipe 


Number of Critical 
Element Velocity 



9.00 


— 2l./5~ 


■3L/5 


a . 65 



L/ 3 L/ 3 -~*»p 


1 L/3 * "' * “ .*! ^ 


11 *76 


,3 ~T— 


-L/2 


9.90 


4 ^ '* ’[a / 3 — Tj/ 2 — *«*~4b ^ 

5 ZT 


11.24 


-%,/6 

6 > KS 



10.26 
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TAP" E 6 

Cr it.ical ve locities i!or Multiple sp a n pip e s w i th one end fixed 

and other supports elas t ic ~5uppo r t s 

P = 0.5, a. = p = p = r = f = O 
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TABLE 7 


Cyltygal velocities for Multi-span verLlc al pipes 
P = 0.5, a « p = p = -r = f = o 



CENTRAf 

Acc.N°- 


URRMW 







TABLE 3 


Values of frequencies bounding p r Inc ipal 
primary instability region for single span 
simply supported pipe 

u = 0.6 it, jS 1 / 2 = 0.8, a =[l = p = y = f = o 
NEL = 5 



• — — — 

Frecruencies 

Lower 

Upper 

0.01 

15.345 

15 .507 

0.10 

14.613 

16.204 

0.20 

13.765 

16 .948 

0 . 30 

12.883 

17.654 

0.40 

11.955 

' 1 

18.319 
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TABLE 9 


frequencies bounding the principal primary 
Instability r egio ns associ a t e d w i th fir st two modes 
of a two span pipe obtained by present analy s is and 

by Kulkarnl f 1 3 ] 


u = 1.2 tt, = o.O, a = P = p = 7 f = 0 


NEL = 8 
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A 



Fig 1 Pip® conveying fluid at velocity V 


A A A A A 



Fig 2 Finite elements of the pipe 





fig. 3 Typical finite element of tfn pipe 










s Regions of principal primary instability of 
one span pipes for 

rl 


IMP l =5 
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Fig, 8 Regions of principal primary instability of a 
two span simply supported pipes for 

Ua0.6'!T,.p t/2 ir 0.4 , oUn*ps?*fsO 8 




Fig, 9 Regions 
span s 


of principal primary instability of th 
imply supported pipes for 




h l / n* 112 H 

a — ~A 

b 2U./4 

.*■" «*““* iv 


Second mod* 

\ _ 




/ 

First mode 

'X- 


Fixed simply 


.*w#C? # \ 


_J_. 

0.1 


Reqions of principal primary instability of two 
span pipes one end fixed and other supports 
simple supports for 
U=1.2TT, p l/Z =0-4,ot- = P = p =■* = f = ° NEL = 




L/34-*L/3 


100 


First mode 



50 L 


"»« SfKWSK ffi^SKSS 

sirnpl« supports for . 

U=I.HTT,p'' 2 =O.A.o(. = r=P = i =‘ i0 NEL = 





I* ‘ ^ ^ 

^r/8.^!.y.8:rL 


Second mode 


First mode 


« h »»« * *mm . ; rzzJt zzz ^" 7 7 ^* 7 ^. ^ 

Fixed -Fixed 


w 

Fig.12 Regions of principal primary Instability for two 
span pipes both ends fixed and other support 
simple support for 

u=1.2TT, p' /2 =0.4,oC=r=p = i = f=o NEL = e 






Second mode 





Fig. 15 Regions of principal Pr'^y ^ability of thr* 
span pipes one end fixed oiner ir«« « 
intermediate supports simple supports or 


Lt-l.alT, P^^«0'4 * oCsPsps^ 


NEL=12 



APPENDIX-1 


The values of various matrices are as follows 


mf 
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-6h 2h 
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3h 
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3h 1 

3h 

v 2 

- h 

36 

- 3h 
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■r x . 
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30h 


36 3h 

- 36 
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36 

6h 

- 36 

0 

4h 2 

- 3h 

.2 

-h 

■r 


2 
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4h 
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6h 2 
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j 





- 


( Al * 3 ) 




69 


APPENDIX- 2 

HOUND ARY CONDITIONS FOR STEADY FLOW 
A 2 . .1 Simple Supports : 

For pipe simply supported at extreme ends , the 
boundary conditions will be 

,2 

v/ a o = 0 (A2.1) 

ar 

For a intermediate simple support, the boundary 
condition will be 


w *= 0 (A2 . 2 ) 

A 2, 2 Fixed Supports : 

For fixed ends of the pipe, the boundary conditions 

will be 

(A2.3) 


M = 0 


w = 0 

A2.3 Plastic Supports , Displacement Springs : 

For pipe supported by a displacement spring, 
3 


JLUJ + a. w = 0 at 5=0 
.3 d 


H' 

a 3 w 

Dfi3 


a d w = 0 


at $ = h 


(A2.4) 


where 


a 


'd 


K cl L ‘ 
E I 


(A2.S) 



70 


A 2. 4 Elastic Supports, Torsional Springs : 

FOr Pip ° flu PPorted by torsional spring, K. v 


o 

ft* V\r 

» ( 2 

v ? 
a w 


a? 


a 


9w 


t = 0 at 5 = o 


71 


+ a 


■ 

fc &£ 


= 0 at K = h 


(A2.6) 


where 


a 


t 


K^L 

eT" 


(A2.7 ) 
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APPENDIX 


I he matrices [ A ] , [b], [e] and [f] can be given as 
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APPENDIX - 4 

HOUND ARY CONDITIONS FOR HARMONICALLY PERTURBED FLOW 
A 4.1 Simple Supports ; 

For pipe simply supported at extreme ends, the 
boundary conditions will be 

Y = 0 2 = 0 

2 

d Y . 2 

— ') ~ 0 d 2 _ n 

d Z A T71 = 0 (A4.1) 

d£ 

For a intermediate simple support, the boundary 
co nd i tic ns wi 1 1 be 

Y = 0 1 Z = 0 ( A 4 . 2 ) 



For a displacement spring K d 


d 3 Y 
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y - o, & + 
a r 
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d 3 Z 

Y * 0, 2-4 - 

dr 

a d z = ° 

at 

tTPf 

a 

(A4.4 ) 
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A 4 . 4 


where 


Elastic Supports, Torsional Springs 
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pipe supported by 

torsional 

spring. 
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